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tasks  in  binary  image  processing.  A  theoretical  model  for 
computation  is  given  using  these  operations  as  basic  instruc¬ 
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within  time  O(logN)  for  input  images  of  size  NxN. 
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1.  Introduction 


In  this  paper,  a  binary  image  is  a  pattern  of  I's  on  a 
background  of  O's,  filling  out  a  scanning  field  of  a  certain 
size  NxM.  In  the  binary  image 

X  =  (X^j)  i=0,l,...,N-l;  j=0,...,M-l 

X^jSO  means  that  (i,j)  is  a  background  point,  and  ^^^=1  that 
(i,j)  is  an  object  point.  It  is  well  known  (see,  e.g.,  [1]) 
that  many  operations  on  binary  images  such  as  edge  detection, 
shrinking,  expanding,  or  projection  can  be  performed  in  parallel. 
The  following  simple  ex^uaple  of  edge  detection  demonstrates  this 
computational  approach: 

Assume  that  all  points  (i,j)  with  X^^=l  represent  a  certain 
object  obtained  as  a  result  of  thresholding  a  gray  level  image 
of  this  object.  The  contour  of  the  object  is  defined  to  be  the 
set  of  all  points  (i,j)  with  X^j=l,  such  that  in  the  neighbor¬ 
hood  of  (i,j)  there  is  at  least  one  point  (k,£)  with  Xj^^=0.  Let 
us  use  the  4 -neighborhood 

ZEZ 

for  points  p  =  (i,j)  in  such  a.  binary  image.  Then,  for  a  binary 
image  X  the  following  procedure  computes  an  edge  image  X*  in 
which  Xjj=l  iff  the  point  (i,j)  is  a  contour  point  in  X  (cp. 
[1,2]) : 

1.  Y  shift  X  one  column  to  right; 

2.  Z  AND(X,Y)  for  all  points  of  X  and  Y  in  parallel; 


3.  Y  shift  X  one  coliinui  to  left; 

4.  Z  ANO(Z,Y)  for  all  points  of  Z  and  Y  in  parallel; 

5.  Y  shift  X  one  row  up; 

6.  Z  ANO(Z,Y)  for  all  points  of  Z  and  Y  in  parallel; 

7.  Y  shift  X  one  row  down; 

8.  Z  AND(Z,Y)  for  all  points  of  Z  and  Y  in  parallel; 

9.  Z  *■  NON(Z)  for  all  points  of  Z  in  parallel; 

10.  X*-*-  AND(X,Z)  for  all  points  of  X  and  Z  in  parallel. 

Thus,  the  computation  of  X*  =  edge (X)  can  be  done  within  10 
steps,  independently  of  the  size  NxM  of  the  binary  image  X, 
where  in  each  step  a  logical  operation  or  a  shift  will  be  per¬ 
formed  on  the  whole  image  field  in  parallel.  On  the  other  hand, 
using  a  sequential  computer,  the  computation  of  X*  requires 
O(NxM)  computational  steps. 

Using  CCD  or  optical  techniques,  such  a  parallel  procedure 
can  be  implemented  in  hardware  [3,4].  However,  our  considera¬ 
tions  are  from  a  theoretical  standpoint  and  are  not  concerned 
with  questions  of  implementation. 

The  paper  is  organized  as  follows:  In  Section  2  the  compu¬ 
tation  model  for  parallel  operations  on  binary  images  is  defined, 
following  the  definition  in  [2].  In  Sections  3,  4,  and  5  we 
discuss  some  examples  of  binary  image  processing  demonstrating 
the  utility  of  the  available  parallelism  for  fast  solutions. 
Section  3  discusses  local  operations;  Section  4  the  recognition 
of  rectangles,  squares,  and  isosceles  right  triangles,  and 


Section  5  the  computation  of  area  and  perimeter  for  a  simply 
connected  object  [1]. 
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2.  Definition  of  the  PBS 


For  reasons  of  simplicity  we  will  consider  bineury  images 

p 

of  size  NaM32  ,  for  integer  p22.  This  choice  is  motivated  by 
the  fact  that  a  power  of  2  is  optimal  for  the  image  size  in 
some  sense:  Binary  image  transforms  via  the  image  processing 
computer-  iiiised  in  this  paper  require  in  many  cases  the  same  time 
qoiplexi^y '.tor  all  sizes  N,  for  any  given  positive 

iinteger  ‘p'.  ' 

^  i  -I  <  *  • 

Let  6  be  the  set  of  all  binary  images  of  size  NxM.  For 
XfcB,  or  X(i,j)  denotes  the  value  of  the  binary  image  X  at 
-the  point  (i,j).  Points  (i,j)  outside  the  scanning  field 
{0,1,.. « .  ,N-1}  are  always  assumed  to  have  value  0. 

Our  model  for  computation  represents  a  hybrid  system  of  a 
vector  machine  [5]  with  a  matrix  machine  [6].  Generally  speak¬ 
ing,  the  matrix  machine  computes  binary  image  transforms,  and 
the  operation  of  this  two-dimensional  machine  is  controlled  by 
switching  in  a  vector  machine.  Both  machines  are  parallel  pro¬ 
cessors  working  as  single- instruction-multiple-data  (SIMD) 
computers.  We  denote  our  model  for  computation  by  PBS  (£aralleles 
Binarbildverarbeitungssystem) 

The  PBS  potentially  makes  use  of  countably  many  registers 
of  three  types; 


^In [2] this  system  was  called  IPC,  but  we  think  that  the  deno¬ 
tation  PBS  is  more  appropriate.  In  fact,  the  definition  of  the 
PBS  fully  agrees  with  the  IPC  definition  given  in  [2] . 


(i)  vector  registers 

(ii)  index  registers 

(iii)  matrix  registers  X,Y,Z,... 

All  registers  may  be  siibscripted. 

Vector  registers  store  one  vector  of  the  set  V;»{0,1}  at 

a  time.  For  simulating  operations  on  nonnegative  integers,  we 

adopt  the  traditional  convention  of  writing  the  sequence  of  bits 

(elements  of  {0,1})  from  right  to  left.  Hence,  for  the  actual 

N 

value  <A>  of  a  vector  register  A  the  relation  0ai<A>*2  -1  is 
always  true.  Furthermore,  we  can  use  vector  registers  for  simu¬ 
lating  operations  on  negative  integers  or  on  real  numbers .  In 
such  cases,  we  need  special  conventions  for  the  meaning  of  some 
bit  positions.  Note  that  the  length  N  of  a  stored  vector  is 
identical  to  the  size  N  of  the  binary  images. 

Index  registers  store  one  vector  of  the  set  0*UO*10*  of 
length  log  N+1  at  a  time.  Hence,  for  the  actual  value  <I>  of 
an  index  register  I  we  have  <I>€{0 , 1, 2, 4 , . . . ,N} ,  at  any  time. 
Usually  these  auxiliary  registers  will  be  used  for  storing  shift 
distances,  or  for  loop  counting. 

Matrix  registers  store  one  binary  image  X€B  at  a  time.  That 

2 

is,  each  register  of  this  type  stores  exactly  N  bits  in  a  two- 
dimensional  square  array.  For  excunple,  in  the  uppermost  row  of 
such  a  register  X  the  N  bits 

^00  ^01  ^02  *  •  *  ^0,N-1 

are  stored  in  this  order. 


The  instructions  for  a  PBS  program  are  drawn  from  the  in¬ 
struction  set  given  in  Table  1.  In  what  follows,  we  will  con¬ 
sider  each  instruction  to  have  a  cost  of  one  unit  of  time 
(uniform  cost  criterion) .  The  time  of  a  computation  is  its 
length  (number  of  executed  basic  instructions) .  For  exaunple,  we 
may  utilize  while- loops,  or  other  bitwise  parallel  Boolean 
operations  such  as  v  and  However,  it  should  be  clear  how 
to  translate  such  constructs  into  the  austere  language  described 
in  Tcdsle  1  at  a  cost  of  a  constant  factor  in  time. 

Note  that  our  instructions  guarantee  that  the  contents  of 
a  vector  register  are  never  assigned  to  an  index  register  and 
vice  versa.  Any  index  register  contains  at  most  one  1  at  any 
time. 

We  write  programs  in  an  Algol-like  notation,  cp.  [7] .  Be¬ 
cause  it  is  not  necessary  that  our  language  should  be  implemented, 
we  take  liberties  with  the  progreunming  language. 

A  program  for  the  PBS  is  a  sequence  of  instructions.  Our 
programs  are  deterministic.  On  the  top  of  a  progreun  we  enumerate 
those  registers  which  will  be  used  during  the  computation.  At 
the  start  of  a  program  all  non-input  registers  are  identically  0 . 

For  exeunple,  the  progreun  for  computing  an  edge  image  (cp. 
the  exaunple  in  the  Introduction)  has  the  form 
begin  index  I;  image  X,Y,Z;  read  X;  int  1=1; 

Y;=X-^I;  Z;=XAY;  Y;=X-^I;  Z;=ZAY; 

Y;=X+I;  Z:=ZAY;  Y:X+I;  Z;=ZAY; 

Z:=^;  Y:=XAZ;  print  Y 


end 


Instruction 


int  A»m;  int  I»2^; 
image  X=Xq; 


read  A;  read  X; 


print  A;  print  X; 


A;=B;  I:«J;  X:=Y; 

A:*B;  A;aBAC; 
X:*Y;  X:=YAZ; 
A:*B-^I;  A;=B-^I; 
I;=J-*-; 

X:=Y-<-I;  X;=Y-^I; 
X:=YfI;  X;=Y4-I; 


if  A{^}0  then. . .else. 
if  l{^}0  then. . .else. 
X:»A; 


A;*X; 


Function 


constant  vectors  or  binary  images 
registers  A,I,X,  esp.  <A>«m,  <I>=2  , 
and  X=Xq68; 

input  statements  for  vector  and 
matrix  registers; 

output  statements  for  vector  and 
matrix  registers; 

direct  assignments; 

bitwise  parallel  Boolean  operations, 

shift  operations,  shift  distances  given 
by  <I>  for  vector  and  matrix  register, 
for  index  registers  the  shift  distance 
is  restricted  to  one  bit  position, 
the  bits  shifted  out  are  discarded, 
the  vacated  positions  are  filled  with 
O's; 

. f i;  test  instructions  for  vector  and 
index  registers; 

.  fi; 

A  into  the  lowest  row  of  X,  the 
vacated  positions  in  X  are  filled 
with  O's; 

the  lowest  row  of  X  into  A. 


Table  1.  Basic  PBS  Instructions. 


f 


According  to  a  given  program  the  PBS  computes  a  certain 
function 

n« 

f;  1/  X  8  2  -^1/  ^  X  8 

Possible  interpretations  of  such  functions  are  listed  in  Table 
2.  Such  a  function  belongs  to  the  complexity  class 
PBS  -  TIME(T(N)) 

iff  this  function  can  be  computed  on  PBS  within  time  0(T(N)). 

In  this,  the  time  complexity  is  taken  as  the  maximum  complexity 
over  all  inputs  (worst-case  complexity) .  E.g.,  we  have  seen 
that  the  function  edge; B-»-B  belongs  to  the  complexity  class 
PBS-TIME(l) . 

According  to  Table  1,  test  instructions  for  matrix  registers 
(Is  X  identical  to  0  in  all  bit  positions,  or  not?)  may  be  per¬ 
formed  by  a  PBS  subroutine  within  time  Odog  N)  .  To  see  this, 
let 

zero (X) ;=if  X=0  then  1  else  0  f i . 

The  function  zero;  belongs  to  the  complexity  class 

PBS-TIME(log  N)  ; 

begin  vector  A;  index  I;  image  X,Z; 
int  I=N/2;  read  X; 
while  1^0  ^ 

Z;=X+I;  X;=XVZ;  od; 


A;=X; 


Function 
8”-^  8 
8"-^  1/ 

(/ 

1/"^  8 


Fields  of  Application 
image  processing 

feature  extraction,  image  coding, 
recognition 

classification,  auxiliary  operations 
image  reconstruction,  image  decoding 


Table  2.  Some  Function  Interpretations. 


The  principle  of  this  prograun  is  easy  to  understand.  This 
example  demonstrates  the  utility  of  power-of-2  shifts  using  a 
divide-and-conquer  approach. 

In  the  next  sections  we  will  consider  some  examples  of 
binary  image  transforms.  These  examples  are  selected  for  pre¬ 
sentation  without  any  special  motives.  We  want  to  emphasize 
the  utility  of  the  available  parallelism  for  image  processing 
tasks.  Some  examples  of  binary  image  processing  using  the  PBS 
can  be  found  in  [2],  e.g.: 

1.  Let  nximber  (X)  be  the  quantity  of  I's  in  the  binary 
image  X.  This  function  nuxnber ;  S->- 1/  belongs  to  the  class  PBS-TIME 

(dog  N)  2)  . 

2.  Let  pro j (X)  be  the  projection  of  X  in  the  column  direc¬ 
tion,  i.e.,  all  I's  in  X  are  shifted  down  to  the  lowest  row. 

This  function  pro  j ;  B-»B  can  be  computed  on  PBS  within  time  0(N). 


3 .  Local  operations 


Let  £:8-^B  be  a  local  operator  on  binary  images,  i.e.,  £ 
is  de£ined  using  a  particular  neighborhood,  and  a  particular 
logical  £unction  on  the  points  in  this  neighborhood.  For 
excunple,  the  operator  edge  is  de£ined  using  the  neighborhood 

^2 

X3  Xq  x^ 


and  the  Boolean  £unction 

XoA(Xj^Vx2Vx3Vx^)  =  XQAnon(Xj^AX2AX3Ax^)  . 

In  general  [2],  an  operator  £:8-^8  is  called  a  local  operator 
o£  degree  n  i££  there  exists  a  totally  ordered  "shi£ting  set" 

E  with  card  E=n,  containing  integer  tuples,  and  there  exists  a 
Boolean  £unction  £*  :  {0 ,  l}”'*-{0 , 1} ,  with  the  £ollowing  property: 

The  set  E  can  be  represented  by  the  n-tuple 
(i2 ,  j  2)  »  •  •  w  (ij^/ jjj)  ]  »  according  to  the  total  order  on  E.  Then, 
via 

£(X)(i,j)  =  if  (i,j)€{0,l,2,...,N-l}^  then 


^i+i2^j+j2 . 

else  0  fi; 

the  operator  £  can  be  computed  by  using  the  Boolean  function 
£*  in  the  environment  E. 

Let  LOC  be  the  class  of  all  local  operators  of  degree  n, 

00 

for  niO,  and  let  L0C:=  U  LOC^  be  the  class  of  all  local  operators. 

n=0 

In  [2]  it  was  shown  that 


1.  if  n<(2N-l)^  LOC^  c  LOC  ; 

n  n+i 

2.  LOC^»LOC,  for  all  na(2N-l)2;  and 

3.  the  cardinality  of  the  class  LOC  is  equal  to 

^2  •  [N*2^-  (N-2”“^-N+l)  +  (N-l) 

From  property  3  we  have  i..  particular  that  there  exist  operators 
f!.8-^8  which  are  not  in  LOC. 

Let  f€LOC  be  an  arbitrary  operator  with  environment 
[  (il»  jl)  /  (13/ j2)  '  •  •  • »  ^  and  with  Boolean  function  f*. 

Then  f  can  be  computed  on  PBS  according  to  the  following  general 
progreun: 

begin  image  X,2j^,2^. . . .  .Z^^;  read  X; 
for  t:=l  step  1  until  n  do 
Z^:»(X-*-j^)  od; 

X:=f*(Z3^,22, 

print  X 

end 


where  for  negative  integers  or  i^  we  use  and 

'M^=+(-i^),  respectively.  This  means  that  if  n  is  a  small 
number,  and  if  the  tuples  (i^,j^)  are  close  to  (0,0),  for  t  * 
l,2,...,n,  then  the  local  operator  f  can  be  computed  on  PBS 
within  constant  time.  The  following  exaunples  of  local  operators 
demonstrate  this: 

(i)  smoothing: 


for  X€B.  Then,  let 


smooth (X) (i, j) and  a(X) (i, j) a2)  or 

(Xj^j*0  and  a(X)(i,j)i3) 
then  1  else  0  li. 

This  local  operation  smooth;  B-*-8  has  the  environment  [(0,0), 
(0,-1) , (0,+l) , (-1,0) , (+1,0) ]  and  the  Boolean  function 

*1  ^*2*5  vx2X2X^XgVx2X2X^Xg  ] 

Vx^Xg (X2*X2) VX2X3 (x^®Xg) vx2X2X^Xg  . 

(ii)  Shrinking;  We  use  the  function  a  defined  adjove. 
Then,  let 

shrink  (X)  (i,  j)  :=if.  and  a(X)  (i,j)»4 

then  1  else  0  f i . 

This  local  operation  shrink;  8->-8  has  the  seune  environment  as 
smooth  and  the  Boolean  function 

X1X2X3X4X5  . 

(iii)  Expanding ;  for 

expand (X) (i, j) ;=if  X^j=l  or  (X^j=0  and  a(X) (i, j) il) 
then  1  else  0  ^ 

we  have  the  same  environment  as  for  smooth  and  the  Boolean 
function 

XiVXi (x2Vx3VX4Vxg) . 

Thus  we  see  that  smooth,  shrink,  expand  are  functions  in  the 
class  PBS-TIME(l). 


I 


4 .  Some  recognition  problems 

In  this  section  we  will  consider  the  recognition  of  rectan¬ 
gles,  squares,  and  isosceles  right  triangles. 

4.1  Rectangles 

Let  RECTcS  be  the  set  of  all  images  which  contain  a  single 
solid  rectangle  of  I's  (with  sides  parallel  to  the  sides  of  the 
image)  on  a  background  of  O's.  Let  rect:  S-*-V  be  the  characteristic 
function  of  RECT.  Because  in  the  images  in  RECT  no  noise  is 
present  the  best  approach  for  computing  rect  seems  to  be  the 
consideration  of  "significant  points"  in  an  input  image  X68. 

This  provides  a  fast  method  for  information  reduction. 

We  will  compute  rect  by  counting  the  local  property  "corner." 
Such  an  approach  can  be  found  in  [8]  using  perceptrons  as  models 
for  computation. 

Proposition  1.  The  set  RECT  can  be  recognized  on  PBS 
within  time  O(log  N) . 

Proof;  We  use  the  patterns 


In  a  first  step  of  the  algorithm,  for  nsl,2,..., 
the  corner  point  p  of  a  binary  image  iff 


8  write  1  at 


(or  a  90**  rotation  of  this)  is  the  place  of  the  observed  pat¬ 
tern  in  the  input  image  X;  at  all  other  points  of  write 
0.  An  example  of  this  operation  is  shown  in  Figure  1.  This 
first  step  is  a  sequence  of  8  local  operators  and  can  be  done 
within  constant  time. 

The  input  image  X  belongs  to  R£CT  iff  each  of  the  images 
^l'^2'^3'^4  exactly  one  point  with  value  1,  and  the  images 
Yg,Yg,Y^,Yg  are  identically  0  at  all  points. 

In  the  second  step  of  the  algorithm,  let  Y;=YgVygVy^vYg. 

If  Y=0  then  go  to  step  3;  otherwise  STOP  with  "X  is  not  in  RECT." 
According  to  our  algorithm  for  zero  (Section  2) ,  this  step  may 
be  performed  within  time  O(log  N) . 

In  the  third  step  of  the  algorithm,  we  check  if  all  of  the 
images  Yj^,Y2/Y2,Y^  have  exactly  one  point  with  value  1.  If  this 
is  true  X  belongs  to  RECT,  otherwise  not. 

The  open  question  for  the  algorithm  just  described  is 
whether  there  exists  a  PBS  program  which  recognizes  if  an  input 
image  X  has  exactly  one  point  with  value  1  within  time  Odog  N)  . 
Let 

one(X) ;=if  X  has  exactly  one  point  with  value  1 
then  1  else  0  fi . 

The  following  program  computes  the  function  one; B-»-l/  within 
Odog  N)  steps: 


input  image  X 


0  0  0  0 
10  10 
1110 
0  10  0 


0  0  0  0 
0  0  0  0 
0  0  10 
0  10  0 


0  0  0  0 
10  10 
0  0  0  0 
0  0  0  0 


0  0  0  0 
10  10 
0  0  0  0 
0  0  0  0 


0  0  0  0 
0  0  0  0 
10  0  0 
0  10  0 


0  0  0  0 
0  10  0 
0  0  0  0 
0  0  0  0 


0  0  0  0 
0  0  0  0 
0  0  0  0 
10  0  0 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  10 


0  0  0  0 
0  10  0 
0  0  0  0 
0  0  0  0 


Figure  1.  Examples  of  images. 


begin  vector  A,B,C;  index  I;  iaege  X,Y,2; 
int  l-N/2;  read  X;  iaiage  Y50; 
while  ^ 

Z;»(X+I)VX;  Y:-YV{(X+I) AX) ; 

X:*Z;  I;»I-*-od; 

A:»X; 

if  A=Q  then  print  0 
else 

if  Y^O  then  print  0 
else  int  I=N/2;  int  C=0; 
while  do 

B;=(A->-I)  VA;  C;=CV(  (A-^1)  AA)  ; 

A:®B;  I;=I-*-od; 
if  Ci^O  then  print  0 
else  print  1  fi; 
fi; 

li? 

end 

After  the  first  while-loop  of  this  program  we  put  the  lowest 
row  of  X  into  the  vector  register  A.  If  in  the  input  image  X 
there  was  any  point  with  value  1  then  A/0/  otherwise  we  stop 
with  output  0  after  O(log  N) steps.  If  in  the  input  image  X  there 
was  any  column  with  more  than  one  point  with  value  1  then  Y/IO, 

The  test  "Y^50?"  needs  O(log  N)  steps.  If  Y^O,  we  stop  with 


I 


output  0  after  O(log  N)  steps.  Now  we  know  that  and  in 
each  column  of  the  input  image  X  there  was  at  most  one  point 
with  value  1.  After  the  second  while- loop  we  check  the  register 
C.  If  in  A  at  the  starting  point  of  the  while- loop  there  was 
more  than  one  bit  position  with  value  1  then  Cj^O ,  and  we  stop 
with  output  0  after  O(log  N)  steps.  If  C=0,  then  in  the  input 
image  X  there  was  exactly  one  column  with  exactly  one  point  with 
value  1.  We  stop  with  output  1  after  O(log  N)  steps.  □ 


4.2  Squares 

Let  SQUAcS  be  the  set  of  all  images  which  contain  a  single 
solid  square  of  I's  (with  sides  parallel  to  the  sides  of  the 
image)  on  a  background  of  O’s.  Let  squat B'^B  be  the  characteristic 
function  of  SQUA. 

Proposition  2.  The  set  SQUA  can  be  recognized  on  PBS  within 
time  O(log  N)  , 

Proof:  According  to  the  relation  SQUACRECT  and  to  Proposi¬ 
tion  1  we  are  able  to  assume  that  the  input  image  X  is  in  RECT 
already.  The  problem  is  to  decide  between  squares  and  nonsquare 
rectangles. 

For  the  computation  of  squa  we  use  the  images  Yj^,Y2»y2'^4 
as  defined  in  Section  4.1.  Let  X;*Yj^vy2VY2VY^;  then  in  X  exactly 
the  four  corners  of  the  input  rectangle  are  the  points  with  value 
1.  In  the  image  Yj^  there  is  a  1  in  the  position  of  the  lower 
right  corner  of  the  rectangle.  We  spread  this  1  along  a  diagonal 
path  toward  the  upper  left  corner  of  the  image.  In  the  image 
Yj  there  is  a  1  in  the  position  of  the  upper  left  corner  of  the 
rectangle.  We  check  if  this  1  is  on  that  diagonal  path  or  not. 

In  the  positive  case  we  know  that  the  rectangle  is  a  square: 
begin  image  Yj^,Y2;  index  I; 
read  Yj^,Y2;  int  1*1; 
while  lj<0  ^ 

Yi:»YiV(  (Yi-^I) +1)  ;  I;«I'^od; 
if  Yj^^O  then  print  1  else  print  0  f i 


end 


This  program  runs  within  time  Odog  N)  .  Note  that  the  1  in 
the  index  register  I  is  shifted  out  of  the  log  N+l  bit  posi¬ 
tions  of  I  to  the  left  after  log  Ntl  runs  through  the  while- 
loop.  □ 

Propositions  1  and  2  show  that  the  sets  RECT  and  SQUA  can 
be  decided  within  the  same  asymptotic  time  as  is  needed  when 
using  UFCA's  as  models  for  computation  [9]. 


4. 3  Isosceles  right  triangles 

Let  IRTRC5  be  the  set  of  all  images  which  contain  a  single 
solid  isosceles  right  triangle  of  I's  (with  two  sides  parallel 
to  the  sides  of  the  image)  on  a  background  of  O's.  Figure  2  _ 
shows  some  ex^unples  of  such  triangles .  Let  irtr ;  B-*-B  be  the 
characteristic  function  of  IRTR. 

Proposition  3.  The  set  IRTR  can  be  recognized  on  PBS  within 
time  O(log  N) . 

Proof:  An  image  X  is  in  IRTR  iff  for  ial,2,3,4,  one  of 
the  Y^-images  marks  the  hypotenuse;  two  of  the  Y^-images  mark 
the  corner  points  of  the  hypotenuse;  one  of  the  Yj^-images  marks 
the  top  opposite  the  hypotenuse;  in  Yj^,Y2/Y2fY^  there  are  no  other 
points  with  value  1;  for  j=s5,6,7,8,  three  of  the  Yj -images  are 
identically  0,  while  the  fourth  Y^-image  marks  the  points  below 
or  above  the  hypotenuse;  and  in  this  fourth  Yj -image  there  is 
no  further  point  with  value  1  (if  the  input  image  is  the  smallest 
triangle  consisting  of  a  single  1  only  the  fourth  Yj-image  is 
identically  0) . 

Figure  3  provides  a  visualization  of  this  property.  Accord¬ 
ing  to  this  property,  the  following  algorithm  is  a  decision  pro¬ 
cedure  for  IRTR: 

Step  1.  Compute  the  images  Yj^,Y2»  . . . » Yg  in  0(1)  steps. 

Step  2.  Check  that  three  images  of  Yg,Yg,Yy,Yg  are  identi¬ 
cally  0;  let  Z  be  the  fourth  one.  This  may  be  performed  within 
time  0(log  N) . 


Ill 
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Figure  2.  Examples  of  images  in  IRTR. 


Input  image  X 


1111 
1110 
110  0 
10  0  0 


0  0  0  1 
0  0  10 
0  10  0 
10  0  0 


0  0  0  1 
0  0  0  0 
0  0  0  0 
0  0  0  0 


10  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


0  0  0  0 
0  0  0  1 
0  0  10 
0  10  0 


Figure  3.  Y^-images  for  an  IRTR  image. 
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step  3.  Check  that  in  three  images  of 
exactly  one  point  with  value  1;  let  these  three 

images.  This  can  be  done  within  time  O(log  N) . 

Step  4.  Check  that  Z  and  have  the  right  connection  to 

each  other,  i.e.,  in  there  are  points  marked  one  diagonal 
below  or  adsove  the  points  which  are  marked  in  Z,  and  in  Y^  there 
is  one  point  with  value  1  less  than  in  Z.  To  do  this,  we  shift 
Y^  one  row  up  or  down,  compare  the  result  with  Z,  using®,  and 
check  if  in  this  resulting  image  there  is  exactly  one  point  with 
value  1. 

S tep  5 .  Check  which  of  Yj^,Y2,Y2  mark  points  in  the  same 
area  as  Y^  (the  corner  points  of  the  hypotenuse)  ;  let  '£^,^2  be 
these  two  images.  This  can  be  done  within  time  O(log  N) ,  using 
®  and  zero. 

Step  6.  Check  if  the  connection  between  the  two  points  in 
Yj^  and  Y2  marks  the  saune  path  as  Y^  within  time  O(log  N)  : 

a)  Spread  the  points  with  value  1  in  Y^^  and  Y2  in  the  dia¬ 
gonal  direction  using  power-of-2  shifts. 

b)  Compare  Y^  with  the  resulting  images  using  A  and  zero. 

Step  7.  Spread  the  points  with  value  1  in  Y^^  and  Y2  in  the 

row  and  column  direction  using  power-of-2  shifts.  Check  that  the 
point  which  is  marked  in  Y^  is  a  crossing  point  of  these  rows 
and  columns,  using  a  and  zero.  This  can  be  done  within  time 
0 (log  N) .  C 


rr  T 


Altogether,  we  have  a  decision  procedure  for  IRTR  running 
within  time  O(log  N) .  Possibly  this  procedure  can  be  simplified 
by  a  particular  property  of  the  Y^-images,  e.g.,  it  may  be  pos¬ 
sible  that  we  have  the  correct  answer  already  after  Step  6. 

Such  an  analysis  to  find  the  fastest  way  (within  the  O(log  N) 
complexity)  would  be  necessary  for  an  implementation  of  this 
procedure . 


5 .  Area  and  perimeter 


Let  S  be  a  simply  connected  object  [1]  in  a  binary  image 
X,  i.e.,  an  object  without  any  holes.  There  exist  several  ways 
to  define  the  area  and  the  perimeter  of  such  an  object.  Conven¬ 
tionally  [1],  the  area  of  S  is  defined  as  the  number  of  points 

in  S.  Using  this  definition,  using  number (X)  we  already  have  the 
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area  of  S  within  O((log  N)  )  steps.  If  we  define  the  perimeter 
as  the  number  of  border  points  (using  the  4-neighborhood)  of  S 
then  we  can  compute  the  perimeter  as  number (edge (X) )  within  time 
0({log  N) too. 

Now,  we  will  define  the  area  and  the  perimeter  of  S  following 
the  definitions  in  [10] .  Let 
area(X) ;=  j  b+i-1 

where  b  is  the  number  of  border  points  of  S  and  i  is  the  number 
of  interior  points  of  S.  (We  assume  here  that  i>0 . )  This 
definition  of  the  area  of  S  follows  the  classical  theorem  due 
to  PICK  (1899)  on  the  area  of  any  simple  polygon  (whose  edges  do 
not  cross  each  other)  whose  vertices  are  lattice  points. 

For  the  computation  of  area,  we  have  b»number (edge (X) )  and 
i=number (edge (X)#X)  after  O((log  N)^)  steps.  For  division  by  2 
we  need  a  special  agreement  for  the  interpretation  of  the  contents 
of  vector  registers.  To  this  end,  we  fix  the  point  position 
between  the  first  and  the  second  bit  positions  from  the  right 
in  such  a  register.  Division  by  2  represents  a  shift  by  one  bit 


position  to  the  right  only.  Using  the  available  Boolean  opera¬ 
tions  and  shifts  on  vector  registers,  addition  and  s\ibtraction 
can  be  performed  within  time  Odog  N)  [5].  Thus  we  have: 

Proposition  4.  The  function  area  can  be  computed  on  PBS 

2 

within  time  O((log  N)  ). 

For  the  definition  of  the  perimeter  of  S,  we  use  the  follow¬ 
ing  NxN  matrices  P  and  Q: 

P(i,j):»if  (i,j)  is  a  border  point 

then  number  of  border  points  in  the  8 -neighborhood  of 
(i,j)  in  the  horizontal  or  vertical  direction 
else  0  fi; 

Q(i, j) :-if  (i,j)  is  a  border  point 

then  number  of  border  points  in  the  8-neighborhood  of 
(i,j)  in  the  diagonal  direction 
else  0  fi; 

for  i,j  =  0,1,..., N-1.  Figure  4  gives  a  visualization  of  these 

two  matrices  P  and  Q.  Then  the  perimeter  of  S  is  defined  as 
,  N-1  N-1  N-1  N-1 

peri(X) ;=  i[  Z  Z  P(i,j)  +  /!•  Z  Z  Q(i,j)]. 
i=0  j=0  i=0  j=0 

For  each  (i,j)  the  values  P(i,j),  Q(i,j)  are  in  the  set  {0,1,2}. 
Therefore,  we  can  represent  P  and  Q  by  two  binary  images 
and  Qj^,Qg,  respectively  (cp.  [6]),  e.g.,  Pj_  (i,  j )  Pg  (i,  j )  denotes 
the  binary  representation  of  P(i,j).  For  the  computation  of 
the  matrices  P  and  Q,  i.e.,  for  the  binary  images  Pj^,Pg,Qj^,Qg 


input  image  X 


00000000 


peri  (X)  s'S+lO/? 


matrix  P 


Z  ZP(i, j)=16=2-8 

i  j 


matrix  Q 


Z  ZQ(i, j)=20=2-10 
i  j 


0  10  1-10  0  0 

•\  /  \ 

0  1111  1^0  0 

0  0^1  1  1  1  1-1 
I  I 

00111111 
'  / 

0  0  1111-10 

0  l^^l  1  /o  0  0 

0  o'^l-l'^o  0  0  0 


0 

0 

0 

0 

0 

0 

0 

0 


0  0  0  0  0  0  0 

1  0^1-1  000 

1  0  0  0  0,0  0 

0  1  0  0  0  1-2 

•  I 

0  2  0  0  0  0  1 

*  / 

0  1  0  0  1-10 

/  / 

0  0  0  0  0  0  0 

\  / 

0  1-1  0  0  0  0 


00000000 

0101-1000 
»\  /  \ 

01200200 
V  \ 

0010001-0 
I  I 

00000001 

'  / 

0  0  1  0  0  1-1  0 

/  / 

02002000 

\  / 

0  0  1-1  0  0  0  0 


Figure  4.  The  matrices  P  and  Q 


we  perform  four  local  operations  on  the  input  image  X  only 
(cp.  [10])  within  constant  time.  Using  the  function  number 
and  the  addition  on  vector  registers  we  get  the  values 
jEZP(i,j)  and  izEQd^j)  within  time  0{(log  N)^).  For  multipli¬ 
cation  by  /7  we  use  the  convention  that  the  first  bit  position 
of  a  vector  register  (i.e.,  on  the  right  end)  indicates  if  the 
number  in  this  vector  register  is  a  multiple  of  /?  or  not.  For 
example 

8 

10^ 

is  the  final  output  for  the  object  S  in  Figure  4.  Thus,  we 
have: 

Proposition  5.  The  function  £eri (X) »a+b/7  can  be  computed 

2 

on  PBS  within  time  0((log  N)  ). 


0|0|0|1|0|0|0|0 


0  0  0  1  0  1  0  1 


6 .  Conclusions 


We  have  defined  a  model  for  parallel  processing  (PBS)  which 
formalizes  Boolean  operations  and  shifts  on  binary  images 
(Boolean  matrices) .  It  was  shown  that  power-of-two  shifts  are 
very  useful  for  fast  programs  using  a  divide-and-conquer  approach. 
The  matrix  registers  of  the  model  used  for  computation  represent 
bit  planes  of  a  possible  hardware  implementation  of  this  system. 
Using  several  bit  planes,  grayscale  pictures  can  be  encoded  in 
this  system.  In  this  sense  the  algorithms  demonstrated  on  binary 
images  (local  operations,  recognition  of  geometric  objects,  area, 
perimeter)  are  only  small  examples  of  problems  which  can  be 
solved  using  the  PBS. 
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